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ABSTRACT 

We use empirical techniques to interpret the near-infrared (near-IR) colours of a sample 
of 5800 galaxies drawn from Sloan Digital Sky Survey (SDSS) main spectroscopic sample 
with YJHK photometry from the UK Infrared Deep Sky Survey (UKIDSS) data release 
one. Our study focuses on the inner 3 arcsec regions of the galaxies sampled by the SDSS fi- 
bre spectra. We study correlations between near-IR colours measured within this aperture and 
physical parameters derived from the spectra. These parameters include specific star forma- 
tion rate (SFR), stellar age, metallicity and dust attenuation. All correlations are analyzed for 
samples of galaxies that are closely matched in redshift, in stellar mass and in concentration 
index. Whereas more strongly star-forming galaxies have bluer optical colours, the opposite 
is true at near-IR wavelengths - galaxies with higher specific star formation rate have redder 
near-IR colours. This result agrees qualitatively with the predictions of models in which Ther- 
mally Pulsing Asymptotic Giant Branch (TP-AGB) stars dominate the H and /-C-band light 
of a galaxy following a burst of star formation. We also find a surprisingly strong correlation 
between the near-IR colours of star-forming galaxies and their dust attenuation as measured 
from the Balmer decrement. Unlike optical colours, however, near-IR colours exhibit very 
little dependence on galaxy inclination. This suggests that the correlation of near-IR colours 
with dust attenuation arises because TP-AGB stars are the main source of dust in the galaxy. 
Finally, we compare the near-IR colours of the galaxies in our sample to the predictions of 
three different stellar population models: the Bruzual & Chariot (2003) model, a preliminary 
version of a new model under development by Chariot & Bruzual, which includes a new 
prescription for AGB star evolution, and the Maraston (2005) model. 

Key words: galaxies: star formation history, near-IR, dust 



1 INTRODUCTION 

Unlike optical colours, the near-IR colours of star-forming galax- 
ies are poorly understood. The area and depth of sky covered by 
near-IR surveys has been steadily increasing in recent years, so it is 
important to understand the nature of the stars that contribute to the 
light emitted in the wavelength range from 1 to 2 /xm. This is nec- 
essary if near-IR magnitudes and colours are to be used as a tool to 
infer physical information about integrated stellar populations. As 
we will see in this paper, the near-IR colours of local spiral galaxies 
turn out to be very sensitive to stars that contribute almost no flux 
at optical wavelengths and that are difficult to model theoretically. 



One of the first observational studies of the integrated 
optical/near-IR colours of star-forming galaxies was carried out 
by Aaronson (1978), who showed that galaxies lie along a well- 
defined morphological sequence in the UVK colour plane in a 
similar manner to the well-known UBV colour- morphology rela- 
tion. Frogel (1985) studied the UBVJHK colours of 19 late-type 
spiral galaxy nuclei and showed that there was no correlation be- 
tween their UBV and JHK colours. This implied that different 
stars contribute to the blue and the infrared light. It was proposed 
that the blue light was dominated by young stars and that the in- 
frared light was dominated by an old and metal-rich stellar popula- 
tion similar to that found in elliptical galaxies. 
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Stellar population synthesis techniques can be used to predict 
which stars contribute most strongly to the near-IR luminosity of 
a galaxy as a function of its age. The results show that the near- 
IR fluxes of galaxies are dominated by cool stars on the red giant 
branch at ages greater than 1-2 Gyr and by stars located on the 
asymptotic giant branch (AGB) at younger ages. In the last phase 
of their evolution, the thermally-pulsing (TP) phase, low metallic- 
ity AGB stars are extremely bright and dominate the K-band light 
of a galaxy even though they are few in number. This evolution- 
ary phase is very difficult to model and can be a source of consid- 
erable discrepancy between different population synthesis models 
(Maraston et al. 2006, Bruzual 2007). 

It has been proposed that the combination of optical and near- 
IR photometry can break the well-known degeneracy between the 
effects of age and metallicity on the stellar populations of galax- 
ies at optical wavelengths. The idea is that the optical colours are 
primarily sensitive to age, because the main contribution to the op- 
tical luminosity of a galaxy comes from upper main sequence and 
turn-off stars. On the other hand the near-IR colours (H and K 
bands) are mainly determined by stars in the AGB phase at low 
metallicity and in the red giant branch stars (RGB) phase at high 
metallicity. Low metallicity carbon-AGB stars are one magnitude 
brighter in the if -band than oxygen-rich M stars and the ratio of 
carbon-rich to oxygen-rich stars decreases sharply with increasing 
metallicity (Lee et al. 2007). Several studies have used the com- 
bination of near-IR and optical broad band colours to estimate the 
ages and metallicities of star-forming galaxies. Generally a max- 
imum likelihood method is used to infer these parameters from 
the nearly orthogonal age-metallicity grids that are predicted using 
stellar population synthesis models. 

De Jong (1996) analysed radial gradients in the B-V and r-K 
colours of 86 face-on spiral galaxies and tested different hypothe- 
ses to explain their origin. His conclusion was that the colour gra- 
dients are best explained by a combined stellar age and metallic- 
ity gradient accross the disk, with the outer regions being younger 
and less metal-rich on average. Similar analyses have been con- 
ducted by Bell & de Jong (2000) and MacArthur et al. (2004) using 
samples of low inclination galaxies that span a range in Hubble 
type. They used r-H/K versus B-r colours to investigate age and 
metallicity gradients and showed that both age and metallicity are 
strongly correlated with local surface brightness. Studying ellip- 
tical and spiral galaxies, Mobasher et al. (1986) have highlighted 
that near-IR colours are particularly sensitive to metallicity. The 
most recent study has been performed by Wu et al. (2(K)5) using a 
sample of 36 nearby early-type galaxies from the Sloan Digital Sky 
Survey. This study found that metallicity variations are responsible 
for the observed colour gradients in these systems. 

These results have recently been challenged by Lee et al. 
(2007), who show that the orthogonal age-metallicity grids com- 
puted using the Bruzual & Chariot (2003) models become much 
more degenerate when more accurate treatments of TP-AGB stars 
and of convective core overshoot are included. However, even after 
these improvements, the models have difficulty in matching data 
on star clusters of known age and metallicity. Geller et al. (2006) 
also encountered difficulties when attempting to model the near-IR 
colours of interacting galaxies in close pairs. They found a subset 
of galaxies with much redder H — K colours than could be ex- 
plained by simple models and interpreted their results as evidence 
for extremely hot (600-1000 K) dust within compact star-forming 
regions. 

In this paper, we adopt an empirical approach. We combine 
broadband YJHK photometry from the UK Infrared Deep Sky 



Survey (UKIDSS) data release one with optical photometry from 
Sloan Digital Sky Survey (SDSS) fourth data release. We calculate 
optical and near-IR colours within the 3 arcsecond diameter fibre 
aperture sampled by the SDSS spectra. This allows us to study 
correlations between both optical and near-IR colours and a variety 
of physical parameters that have been derived directly from the 
spectra; these include specific star formation rate, metallicity, 
dust attenuation, and mean stellar age. Because UKIDSS and 
SDSS spectroscopy are available for a sample of more than 5000 
galaxies, we are able to study these correlations for subsamples of 
galaxies that are closely matched in stellar mass and in redshift. 
We also compare the optical and near-IR colours of the galaxies 
in our sample to the predictions of the stellar synthesis population 
models of Bruzual & Chariot (2003), a preliminary version of an 
improved model under development by Chariot & Bruzual (2007) 
and the models of Maraston (2005), which implement the TP-AGB 
phase using "fuel consumption" approch. 

This paper is organised as follows. We describe our optical 
and near-IR data in section 2, as well as the physical parameters 
used for our analysis. Our methods for analyzing the correlations 
between near-IR colours and physical parameters are presented in 
section 3 and 4. In section 5, we compare the near-IR colours of the 
galaxies in our sample with the colours predicted by stellar popu- 
lation models. Finally we discuss the implications of our results in 
section 6 and we summarize our findings in section 7. 



2 THE DATA 

2.1 Optical photometry 

The optical photometry is drawn from the Sloan Digital Sky Sur- 
vey (SDSS) photometric galaxy catalogue (York et al. 2000). 
We restrict our sample to galaxies in the publically available 
MPA SDSS database (http://www.mpa-garching.mpg.de/SDSS), 
for which physical parameters have been derived from the galaxy 
spectra. Our optical sample contains 473,034 galaxies from the 
main spectroscopic sample of the SDSS fourth data release (DR4; 
Adelman-McCarthy et al. 2006), covering 4783 square degrees to 
a depth of 17.77 in r-band Petrosian AB magnitude. The magni- 
tudes of interest are measured through the fibre aperture in the u, 
g, r, i, z filters. We have applied corrections for foreground Galac- 
tic extinction according to Schlegel et al. (1998) and we have also 
/^-corrected the magnitudes to redshift 0. 1 (the mean redshift of the 
sample) using v3 of the KCORRECT code (Blanton et al. 2003) 
and assuming a standard cosmology with fi = 0.3, A = 0.7 and 
Ho = 70 km s~ Mpc~ ^ . Our sample covers a redshift range from 
0.005 to 0.3. In the following sections, we will use the concentra- 
tion parameter C, defined as the ratio of the radii enclosing 90 and 
50 per cent of the galaxy light in the r-band, as an indicator of 
the morphological type of the galaxy. The axis ratio of the galaxy 
is defined as the ratio of the isophotal minor axis to the isophotal 
major axis in the r-band, where the isophotes are measured at 25 
magnitudes per square arcsecond. 

2.2 Near infrared photometry 

The near-IR data come from the UK Infrared Deep Sky Survey 
(Lawrence et al. 2007). The Large Area Survey, the largest of the 
five surveys that constitute UKIDSS, is the near-IR counterpart of 
SDSS. The images are obtained by the Wide-Field Camera on the 
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2.3.1 Star-forming and non star-forming classes 

In order to select a sample of star-forming galaxies free of AGN 
contamination, we use the emission-line classifications of Brinch- 
mann et al. (2004). The classifications are based on the location 
of galaxies in the Baldwin, Phillips & Terlevich diagram (Baldwin 
et al. 1981). Galaxies are divided into 6 classes: star-forming, com- 
posite, AGN, low S/N AGN, low S/N star-forming, and unclassifi- 
able. We selected the star-forming and the non star-forming class 
(unclassifiable) for our analysis. The average values of the physical 
parameters are listed in Table 1 for these two classes. As noted by 
Brinchmann et al. (2004), it is important to remember that this is a 
nuclear classification for the galaxies lying at redshift 0. 1 and be- 
low, as these measurements are restricted to the central regions of 
the galaxies. 



Figure 1. Response functions of the UKIDSS YJHK filters (black curves) 
are plotted on top of spectra computed for a star-forming galaxy (blue) and 
for a non star-forming galaxy (orange) at 2: = 0. 1 using the Bruzual & Char- 
lot (2003) models. 



UK Infrared Telescope in Hawaii. At the end of the 7 year observ- 
ing campaign, it will cover 4000 square degrees of the Northern 
sky in the YJHK filters and it will reach a limiting if-band mag- 
nitude that is three magnitudes deeper than that of 2MASS (Skrut- 
skie et al. 2006). In the present study we focus on galaxies con- 
tained in the first data release of the Large Area Survey (Warren 
et al. 2007). This reaches a depth of 18 in if-band Vega magni- 
tude and covers 320 square degrees of the Northern sky. Our near- 
IR sample comprises 2.5 million objects with measured A'-band 
fluxes. The four filters cover the near infrared wavelenths from 
1 micron for the Y filter to 2.4 micron for the K filter. Circular 
aperture photometry is available for 13 diameters. We extract the 
Vega magnitude measured through a 2.8 arcsec diameter aperture 
^, which is well-matched to the 3 arcsec diameter fibre aperture of 
SDSS. These magnitudes are corrected for Galactic extinction and 
fe-corrected to redshift 0.1. In Fig. 1, we show the response func- 
tion of the UKIDSS YJHK filters on top of spectra computed for 
a 10 Gyr old star-forming galaxy (blue) and for a 10 Gyr old non 
star-forming galaxy (orange) at z = 0. 1 using the Bruzual & Chariot 
(2003) models. 

2.3 Physical parameters 

In this section, we briefly define the physical parameters that are 
used in our analysis: specific star formation rate, age, stellar and 
interstellar metallicity and dust attenuation. These quantities are all 
available in the public SDSS MPA database^ and have been derived 
from the galaxy spectra. We note that the fibres in the SDSS spec- 
trograph have diameters of 3 arcsec. At the mean redshift of our 
sample this corresponds to a diameter of 7 kpc (h = 0.7). Analyzing 
galaxy properties within a fixed fibre aperture means that the phys- 
ical size of the region for which we carry out our analysis will be 
larger at higher redshifts. 

^ From the UKIDSS catalogue apermag4 magnitudes were selected, which 
have been corrected for the effects of aperture and seeing under the assump- 
tion that the object is a point source. However, we have verified that the 
resulting colours are unaffected by this correction 
^ http://www.mpa-garching.mpg.de/SDSS/ 



2.3.2 Star formation rate 

The star formation rate (SFR) measured inside the fibre is also 
taken from the Brinchmann et al. (2004) analysis. To estimate this 
parameter, these authors matched the intensity of a variety of emis- 
sion lines, including Ha and H/J, to a grid of HII region photo- 
ionization models. Because the stars that contribute to the ionizing 
flux are massive and very young, with a lifetime < 20 Myr, this 
method provides a measure of the instantaneous star formation rate, 
independent of the previous star formation history. For our study, 
we use the specific SFR measured inside the fibre (SFR/M*). 



2.3.3 Age and metallicity 

We use the light-weighted stellar ages and the stellar metallicities 
estimated by Gallazzi et al. (2005) using a carefully chosen set of 
absorption features that can break the age-metallicity degeneracy 
and that depend weakly on the a/Fe abundance ratio: D4000, H/3, 
H5a+H7a, [Mg2Fe] and [MgFe]'. While the composite Mg-l-Fe 
indices are a good indicator of metallicity, the H-Balmer lines are 
mostly sensitive to age. For star-forming galaxies, we use the gas- 
phase oxygen abundance derived by Tremonti et al. (2004) from 
the strong optical nebular emission lines as the metallicity indi- 
cator, because for these galaxies, the gas-phase metallicity can be 
measured with better precision than the stellar metallicity. These 
values are expressed in 12 -1- log (O/H), where the solar metallicity 
in these units is 8.69. 



2.3.4 Dust 

As an estimate of the dust content of our star-forming galaxies, 
we have used the ratio of the Ha to the H/3 emission line fluxes. 
These emission lines are measured after the subtraction of the con- 
tinuum spectra. As Ha and H/3 are attenuated differently because of 
their different wavelengths, their ratio is a measure of the amount 
of dust present in and around HII regions, where the star forma- 
tion occurs. We assume a value of 2.86 corresponds to the dust-free 
case. An independent measurement of the dust content is provided 
by the z-band attenuation parameter Az derived by Kauffmann et 
al. (2003a), which measures the amount of dust in the interstellar 
medium of the galaxy. It is calculated by fitting spectra with stel- 
lar population synthesis models with an attenuation curve of the 
form ta oc A""'^. The parameter A^ can be estimated for both star- 
forming and non star-forming galaxies. 
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Table 1. Values of the physical parameters for the star-forming and non star-forming samples 





Mfib 


SFR/M* 


Age 


Zgas 


^stellar 


Ha/H/3 






[log M0] 


[logyr-i] 


[Gyr] 


[12+log(0/H)] 


[2©] 




[mag] 


SF: Mean value 


9.6 


-9.9 


2 


9.0 


0.7 


4.1 


0.6 


SF: 1 sigma 


0.7 


0.4 


1 


0.2 


0.5 


0.7 


0.3 


NSF: Mean value 


10.1 




7 




1.1 




0.0 


NSF: 1 sigma 


0.5 




2 




0.4 




0.2 



2.4 Matching near-IR and optical 

We have matched the SDSS and UKIDSS samples in right ascen- 
sion and declination with a maximal allowed separation of 0.5 arc- 
sec in both coordinates. This upper limit is consistent with the mea- 
surement errors on the position of SDSS and UKIDSS, where 96% 
of SDSS objects have an error < 0.2 arcsec in right ascension and 
decUnation and all UKIDSS objects have an error < 0.26 arcsec. 
The mean separation obtained for the match is 0.17 arcsec. The 
matched sample consists of 12,415 galaxies with detections in all 
four near-IR bands. Nearly all SDSS galaxies with spectroscopy 
that are in the field observed by UKIDSS have a near-IR counter- 
part. The star-forming and non star-forming galaxy classes consist 
of 2377 and 3438 galaxies respectively. 



define its partner as the galaxy minimising the quantity sep defined 
by equation (1) below and satisfying the limits in AlogMgb, Az 
and AC defined above, where the denominators are the average 
of the absolute values. After this operation, the two galaxies are 
removed from the list, and we repeat the procedure with the rest of 
the sample. 

^ lAZogMfibl |A^| |Ag| 

P avg(IA«o(;A/Hbl) avgdA^I) ^ avg(|AC|) ^ ^ 

This results in 708 and 673 galaxy pairs for the star-forming and 
non star-forming classes, respectively. The pairs are required to 
have unproblcmatic measurements of mass, specific SFR (for the 
star-forming class only), age, metallicity and dust attenuation for 
both galaxies. 



3 METHODS 



3.2 Calculation of the correlations 



In this section we present our method to derive the correlation be- 
tween the optical/near-IR colours and the physical parameters de- 
scribed in section 2.3. We analyse the star-forming and the non star- 
forming galaxies separately. The pair matching method, described 
in the following section, allows us to suppress the dominant corre- 
lation with stellar mass. For each class, we then evaluate the corre- 
lation between the optical/near-IR colours and specific SFR, stellar 
age, metalUcity and dust attenuation of paired galaxies using the 
Pearson correlation coeficient. 



3.1 Matched galaxy pairs 

It is well known that colours are strongly correlated with galaxy 
mass; the most massive galaxies have redder optical colours and 
the same is true for near-IR colours. The main physical driver of 
this trend is the fact that low-mass galaxies are currently forming 
stars and high-mass galaxies have very little ongoing star formation 
(Kauffmaim et al. 2003b). If we wish to isolate the effects of other 
physical parameters on colours, it is useful to remove this domi- 
nant trend by comparing galaxies with a fixed stellar mass. We use 
a pair-matching technique; for each galaxy, we find a "partner" that 
is closely matched in stellar mass, in redshift and in concentration 
index. We require a close match in redshift, because for a given stel- 
lar mass, the sample is biased towards bluer galaxies at higher red- 
shifts. Matching in redshift also ensures that we compare colours 
within the same physical radius. Matching in concentration means 
that we compare galaxies with the same morphology, because the 
colour is influenced by the presence of a bulge component in the 
galaxy. To summarize, our matched galaxy pairs have: AlogMflb 
< 0.05 with Mfib expressed in M©, Az < 0.02 and AC < 0.1. 

We proceed as follows: we begin by randomising the order 
of the galaxies in our sample. Then we take the first galaxy and 



We use our pair sample to study correlations between the colour 
differences of paired galaxies and differences in specific SFR, stel- 
lar age, metallicity and dust attenuation. The degree of correlation 
is evaluated using the Pearson's product-moment correlation coef- 
ficient r, defined as 



r{x,y) 



(:ov(j:, y) 
a{x) X a{y) 



(2) 



where a is the standard deviation. The correlation coefficient mea- 
sures the degree of linear dependence between the two variables. 
When they are independent, then this coefficient is zero. In gen- 
eral, r varies between — I and -i-l, where r = — 1 means that x and 
y are linearly dependent with a negative slope, and r = +1 means 
that they are linearly dependent with a positive slope. The errors 
on the correlation coefficients are calculated analytically by com- 
bining the measurement errors of the colours and the errors of the 
derived physical parameters (see Brinchmaim et al 2004 for more 
details). We also use a resampling technique to assess whether the 
derived value of r depends on the choice of galaxy pairs. The final 
values of r and their associated errors represent the average of 1000 
resampUngs of the pair-matching process. 



4 RESULTS FOR THE CORRELATION 

4.1 Correlations for star-forming galaxies 

The results of our analysis for the star-forming sample are sum- 
marised in Table2. Given the null hypothesis that they are uncor- 
rected and taking into account the size of the samples, the prob- 
ability that |7-| should be larger than 0.074 is 5% i.e. correlation 
coefficients below this value should be regarded as insignificant. 
The correlation coefficients have been calculated for optical and 



© 0000 RAS, MNRAS 000, 000-000 



Physical interpretation of the near-IR colours of low redshift galaxies 5 



Li_ 



I 

\ 

o 



+ 

CM 



"SI 

T 



8 

I 




-0.4 -0.2 0.0 0.2 0.4 -0.2 0.0 0.2 
A(g - r) [AB] A(H - K) [Vega] 



Figure 2. Correlation between the differences in tlie specific SFR, metal- 
licities, Balmer decrements and axis ratios of tlie galaxy pair and the differ- 
ences in their g — r and H — K colours. The contours indicate the density 
of points. Results are shown for the star-forming sample. 



near-IR colours as a function of specific SFR, mean stellar age, gas- 
phase metallicity, Ha/H/3 ratio, z-band attenuation, and isophotal 
axis ratio. We note that all these parameters except the last one are 
derived from the galaxy spectra. Results are listed for the optical 
colours g — r and i — z and the infrared colours Y — J and H — K. 
The chosen set of colours enables us to study how the shape of a 
galaxy's spectral energy distribution depends on physical parame- 
ters from the optical to the near-IR in a continuous way. We also 
show examples of these correlations in Fig.2: g — r and H — K 
colour differences as a function of differences in specific SFR, gas- 
phase metallicity, Ha/H/S ratio and isophotal axis ratio. 

As can be seen from Table2 and from Figure2, the g~r colour 
is strongly anti-correlated with the specific SFR (r = —0.49), with 
the bluest colours corresponding to the most strongly star-forming 



galaxies. This result is expected, because the emission from young 
stars peaks at UV and blue wavelengths. This anti-correlation still 
exists for the i — z colour, but the H — K colour displays a posi- 
tive correlation with specific SFR (r = 0.23 as seen in Fig.2). At a 
constant mass, strongly star-forming galaxies have redder near-IR 
colours. Even if this correlation is not as strong as the one in the 
optical, it is significant. One possibility is that that AGB stars are 
responsible for this correlation, because they dominate the A'-band 
luminosity of a galaxy 0.1 Gyr after a burst of star formation (see 
section 5.2 for detailed predictions from stellar population mod- 
els). Galaxies with ongoing star formation are rich in AGB stars 
and therefore have redder near-IR colours than quiescient galaxies 
of the same mass. 

The g — r colour also correlates with the mean stellar age mea- 
sured from stellar-absorption lines; galaxies with older ages have 
redder optical colours. This correlation becomes weaker at longer 
wavelengths, and once again the sign of the correlation reverses for 
the H — K colour (r = —0.17). The gas-phase metallicity does not 
seem to correlate with any colour. At a first glance, this is a sur- 
prising result, because metal-rich stellar populations are predicted 
to have redder colours. The most natural explanation is that the 
metallicity of a galaxy is very strongly correlated with its mass and 
not with any other parameter (Tremonti et al. 2004). Because our 
galaxy pairs are constrained to have the same stellar mass, the cor- 
relation with metallicity is strongly suppressed. 

The correlation coefficients for the Ha/H/3 ratio and the z- 
band attenuation are also small for optical colours. Curiously, the 
strength of the correlation increases towards near-IR wavelengths 
(r = 0.25 and 0.23 for the correlation between H — K colour and 
Hq/H/9 and respectively). This result is very surprising, as dust 
should affect optical colours more than near-IR colours. 

One way to gain further insight into this result is to look at the 
correlation between optical and near-IR colours and the axis ratios 
of the galaxies (see Table2 and the bottom panels in Fig.2). Our re- 
sults show that low-inclination (&/a ~ 1) galaxies have bluer optical 
colours than high-inclination {bla <^ 1) galaxies. This is most likely 
due to dust reddening, because the light coming from the center of 
high-inclination galaxies passes through a large amount of mate- 
rial, whereas for low-inclination galaxies it comes directly to us. 
As seen from Table2 the correlation of the near-IR colours with the 
axis ratios is much weaker. This leads us to suggest that the positive 
correlation of the near-IR colours with dust attenuation is caused by 
some phenomenon other than dust reddening. AGB stars are known 
to release large amounts of dust. We thus hypothesize that galaxies 
with redder H — K colours contain a higher fraction of AGB stars, 
and by extension a larger amount of dust. The apparent tendency 
for the near-IR colours to correlate with Balmer decrement and the 
z-band attenuation is thus an induced correlation, not a primary 
one. 



4.2 Correlations for non star-forming galaxies 

Correlations calculated for the non star-forming galaxy sample are 
presented in Table3. We show the correlation coefficient and its er- 
ror for the relations between optical/near-IR colours and mean stel- 
lar age, stellar metallicity, z-band dust attenuation, and the isopho- 
tal axis ratio. All the correlation coefficients are smaller than those 
obtained for the star-forming sample, but the overall results are 
quite similar. The strongest correlation is that between age and 
g — r colour (r = 0.4), and the strength of the correlation de- 
creases at longer wavelengths (r = —0.04 for the H — K colour). 
Table3 shows that the correlation with stellar metallicity is below 
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Table 2. Correlation coefficients for star-forming galaxies for the relations between optical/near-IR colours and the following physical parameters: SFR/M* 
(the specific SFR measured inside the fibre), mean stellar age, Zgas (the gas-phase metallicity), the Ha/H/3 ratio (Bahner decrement), Az (the attenuation of 
the z-band Ught by dust), and bla (the isophotal axis ratio). The errors are estimated as described in the text. 



Correlation Coefficient 


Colour 


SFR/M* 


Age 




Ha/H/3 


Az 


bla 


9- 


r 


-0.49 ± 0.05 


0.62 ± 0.09 


0.09 ± 0.05 


0.11 ±0.03 


-0.09 ± 0.03 


-0.34 ± 0.04 


i — 


z 


-0.37 ± 0.05 


0.46 ± 0.08 


0.06 ± 0.05 


0.13 ± 0.03 


-0.04 ± 0.04 


-0.31 ±0.04 


Y - 


J 


-0.02 ± 0.02 


0.07 ± 0.05 


0.08 ± 0.03 


0.20 ± 0.03 


0.12 ±0.03 


-0.20 ± 0.02 


H - 


K 


0.23 ± 0.03 


-0.17 ±0.05 


0.03 ± 0.03 


0.25 ± 0.03 


0.23 ± 0.04 


-0.12 ± 0.03 



Table 3. Correlation coefficients for non-star-forming galaxies for the relations between optical/near-IR colours and the following physical parameters: mean 
stellar age, stellar metallicity, Az (the attenuation of the z-band Ught by dust), and bla (the isophotal axis ratio). 



Correlation Coefficient 



Colour Age ^stellar Az bla 

g-r 0.4 ±0.1 0.05 ± 0.08 -0.02 ± 0.05 -0.08 ±0.11 

i-z 0.20 ±0.09 0.00 ± 0.07 0.16 ± 0.05 -0.08 ± 0.05 

y - J 0.08 ± 0.04 0.04 ± 0.04 0.02 ± 0.02 -0.04 ± 0.02 

H-K -0.04 ±0.04 0.07 ±0.04 0.14 ±0.03 -0.01 ± 0.03 



the 5% significance level (r ^ 0.07). Once again, the correlation 
of colours with z-hmd dust attenuation increases at near-IR wave- 
lenghs, while the absolute correlation with isophotal axis ratio de- 
creases. We note that the correlation between g — r colour and 
isophotal axis ratio is less significant for this sample, as most of 
the galaxies are ellipticals. 



5 COMPARISON WITH STELLAR POPULATION 
MODELS 

5.1 Description of Bruzual & Chariot (2003) and Chariot & 
Bruzual (2007) models 

We use the Bruzual & Chariot (2003) (hereafter BC03) stellar pop- 
ulation synthesis models to compute the optical and noar-IR colours 
of model "galaxies" for comparison with the data. We also test here 
a preUminary version of the Chariot & Bruzual (2007) (hereafter 
CB07) code, which includes a new prescription for AGB star evo- 
lution. We adopt the Padova 1994 single stellar population evolu- 
tionary tracks and the Chabrier initial mass I'Linclion. The models 
allow us to compute the evolution of the spectrum and the colours 
of a composite stellar population between and 20 G3T for a given 
metallicity and star formation history. For each model galaxy, the 
output colours are evolution and fc-corrected to 2: = 0.1, so that they 
can be compared directly with the observational data. 

We have computed a grid of models with different star for- 
mation histories and metalUcities that are designed to span the ob- 
served colours of the galaxies in our sample. We parametrize an 



' The BC03 model package allows one to choose between Padova 1994 
and Padova 2000 evolutionary tracks, but the latter tends to produce worse 
agreement with observed galaxy colours. On the other hand, as shown by 
Bruzual & Chariot (2003), the choice of the initial mass function should not 
influence significantly the output colours. 



exponential declining SFR as 

0exp(t) = , (3) 

T 

where t is the age and t is the star formation timescale. To repre- 
sent the star-forming sample, we have chosen two models. The first 
one has a constant SFR, and the second one has an exponential de- 
clining SFR with r = 3 Gyr. In general, the star formation history 
of early-type spiral galaxies can be described by an exponential 
law, because the gas available for forming new stars decreases with 
time (Kennicutt 1998). The constant SFR model describes galaxies 
where the available gas is replenished through further infall, and 
is therefore suited to describe later-type galaxies with higher rates 
of ongoing star formation. To model our non star-forming sample, 
we choose r = 1 Gyr. In this model, the star formation rate de- 
creases rapidly and reaches values close to zero a.lt = 5 Gyr. This 
model is designed to reproduce the star formation histories of non- 
interacting elliptical galaxies. 

We generate these models for five stellar metallicities: 

0. 02Zq, 0.2Zq, OAZq, Zq, 2Zq. This range is broader than that 
spanned by the galaxies in our samples, which vary between O.IZq 
and ^Zq. Note that our model galaxies have a fixed metallicity, 

1. e. there is no metal enrichment with age. We plot galaxy colours 
at model ages of 5, 10 and 15 Gyr (we expect the typical age of 
the galaxies in our sample to be around 10 Gyr). The upper limit 
of 15 Gyr is older than currently favored estimates of the age of 
the Universe, but in a few cases models with large ages are needed 
to fit the data. We note that by "age", we mean here the time for 
which the model galaxy has been forming its stars; this must be 
distinguished from the observed age of the stellar population (such 
as that considered in section2.3.3), which is usually a measure of 
its luminosity-weighted mean stellar age. The luminosity-weighted 
age is generally considerably lower, because young stars always 
dominate the optical luminosity in a composite stellar population. 

The effect of dust attenuation on the colours is included using 
the simple two-component model of Chariot & Fall (2000). This ac- 
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counts for different attenuations affecting young stars in tlieir birth 
clouds (i.e., giant molecular clouds) and older stars in tlie diffuse 
ambient (i.e. diffuse) interstellar medium of a galaxy. The main two 
adjustable parameters of this model are the total effective V-band 
optical depth affecting stars younger than 10^ yr (the typical life- 
time of a giant molecular cloud), tv, and the fraction of that optical 
depth that is contributed by dust in the ambient interstellar medium, 
/i. Our models include dust attenuation with ry = 1 and /i = 0.3. 
This corresponds to an attenuation of 1 . 1 magnitude in the V-band 
for stars younger than lO'^ yr and 0.3 mag for older stars. We note 
that the estimated z-band attenuation varies between 0.4 magnitude 
and 0.8 magnitude for the star-forming galaxies in our sample and 
between 0. and 0.03 for the non star-forming galaxies. The atten- 
uation curve across the spectrum has the form t\ = ry (A/5500 
jjjg effect of dust on the colours is indicated by an orange 
arrow on the figures in section 5.3. The length of the arrow corre- 
sponds to a galaxy with f v = 2 and fi = 0.5, corresponding to an 
attenuation of 2.2 mag in the V-band for stars younger than 10*^ yr 
and of 1.1 mag for older stars. This is a reasonable upper limit to 
attenuation by dust. The cross on the arrow indicates the amount of 
reddening that is already included in the models. 
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Figure 3. Evolution of the near-IR colours predicted by the BC03 and CB07 
stellar population models following a burst of star formation. Blue, green, 
yellow, orange and red lines correspond to models with 0.02^0, 0.2^0, 
OAZq , Zq and 2^0 respectively. 



5.2 Difference between the BC03 and preliminary CB07 
models 

The preliminary version of the CB07 model used here differs from 
the BC03 models only in the treatment of AGB stars. The new 
models include the prescription of Marigo & Girardi (2007) for the 
TP-AGB evolution of low- and intermediate-mass stars. This pre- 
scription includes several important theoretical improvements over 
previous calculations, and it has been calibrated using carbon star 
luminosity functions in the Magellanic Clouds and TP-AGB life- 
times (star counts) in Magellanic Cloud clusters (we refer to the pa- 
per of Marigo & Girardi (2007) for detail). As outlined by Bruzual 
(2007), this implies brighter if -band magnitudes and redder near- 
IR colours for the preliminary CB07 model than for the BC03 
model. We caution the reader that the spectra of TP-AGB stars and 
the stellar evolution prescription for pre- AGB phases have not been 
updated in the preliminary version of CB07 model used here. The 
final version of this model will include new spectra for TP-AGB 
stars and an updated library of Marigo & Girardi (2007) TP-AGB 
calculations connected to a new comprehensive grid of pre-AGB 
stellar evolutionary tracks by Bertelli et al. (2007, in preparation). 
Nevertheless we feel that this preliminary analysis is informative. 

To illustrate the differences between the BC03 and CB07 mod- 
els, we have computed the evolution of the near-IR colours after a 
burst of star formation. We assume that one third of the final mass 
of the galaxy is formed during the burst, which occurs at t = 9 Gyr 
on top of a galaxy with exponentially declining SFR. As seen from 
Figure 3, when TP-AGB stars appear, about 0. 1 Gyr after the burst, 
the CB07 models are significantly redder (from a tenth to a few 
tenths of a magnitude) m J ~ H and H — K than the BC03 model. 
This effect is seen mainly at metallicities of 0.2Zq and below. At 
such metallicities, the TP-AGB (and in particular, luminous carbon 
stars) dominate the near-IR light of intermediate-age stellar popula- 
tions in the CB07 model, while RGB stars (which are less luminous 
than TP-AGB stars) prevail at medium to high metallicities. 

5.3 Results from BC03 and preliminary CB07 models 

Figures 4, 5 and 6, 7 present comparisons of the BC03 and pre- 
liminary CB07 models with our data. The observational data are 



plotted as black dots in four colour-colour diagrams: g — r versus 
u ~ g, J ~ H versus Y — J, Y — K versus H — K and g — r 
versus Y — K. The model grids are defined by 3 different ages 
and 5 metallicities from 0.02^0 to 2^0. The dotted, continuous 
and dashed-dotted lines correspond to model ages of 5 , 10 and 
15 Gyr, respectively. For the star-forming sample, the models with 
constant SFR are plotted as stars, while the models with star for- 
mation timescale r = 3 Gyr are plotted as filled squares (Figs. 4 and 
5). For the non star-forming sample, all models have r = 1 Gyr and 
are plotted as filled triangles (Figs. 6 and 7). All models in Figs.4— 7 
include dust attenuation with tv = i and ji = 03. The orange ar- 
row on the lower right comer of each panel indicates the amount of 
reddening for a galaxy with fv = 2 and p = 0.5. The orange cross 
indicates the reddening that is already included in our models. 

5.3.1 Star-forming galaxies 

As seen from Fig. 4, the BC03 models cover the data in the optical 
g ~ r versus u ~ g colour-colour plane. The very bluest galaxies 
in u — g can be fitted if a burst of recent star formation is super- 
posed on the star formation history used in our models. The colours 
of the reddest galaxies can be understood if they are significantly 
reddened by dust. We therefore conclude that our choices of star 
formation history, age and metallicity are appropriate for the star- 
forming sample. These models should then be able to reproduce the 
near-IR colours of the galaxies in our sample. 

We find that, however, the BC03 models fail to cover the data 
in the J — H versus Y — J and Y — K versus H — K colour- 
colour planes. In the near-IR, these models are sensitive primarily 
to metallicity. As can be seen, the models cannot reproduce the 
observed scatter of the star-forming galaxies in the near-IR colour- 
colour diagrams. Neither reddening nor a scatter in age or star for- 
mation history can account for the observed spread. We also note 
that the J ~ H colours of the model tracks are too blue by ap- 
proximatively 0.1 magnitude as compared to the data. This may 
not seem like a large offset, but as the near-IR colours of galaxies 
cover a small range of magnitude (0.4 magnitude in J — H) this is 
not a negligible discrepancy. 

Fig. 5 shows the comparison with the preliminary version 
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Figure 4. Colour-colour diagrams showing BC03 model tracks superimposed on our sample of star-forming galaxies (black dots). The stars and the squares 
denote models with constant star-formation and exponentially declining star formation with r = 3 Gyr. The grids are drawn for 5 metallcities and three ages: 
the coloured lines represent isometalhcity colours of 0.2Zq, 0.4Zq,Zq, 2Zq and the dotted, continuous and dashed-dotted line show isochrones of 5, 10 and 
15 Gyr respectively. The orange arrow indicates the amount of reddening for a galaxy with ry = 2 and /i = 0.5. The orange cross shows the reddening that is 
already included in our models. 
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of the CB07 models. As expected, the resuhs do not differ from 
the BC03 models at optical wavelengths. On the other hand, the 
near-IR colours are now dependent on stellar age at metalUcities 
< OAZq. Although this could in principle provide a mechanism 
for understanding the observed scatter in the near-IR colour-colour 
planes, wc find that these more recent models still cannot fit the 
detailed trends seen in the data. The H — K colours are skewed to 
significantly redder values in the CB07 models as compared to the 
data and the models fail to reproduce the observed spread m J — H 
colour at F — J colours greater than 0.4. This conclusion is likely to 
change when new spectra are adopted for TP-AGB stars in the final 
models. The Y — J and Y — K colours are better reproduced by 
the models. Both these colours are predicted to be good indicators 
of metallicity. In particular, the Y — K colour has almost no depen- 
dence on age in both the BC03 and the CB07 models. We note that 
there is an offset in the UKIDSS Y-band magnitudes of ~0.09 in 
the sense that they are too faint (RHewett, private corrmiunication). 
Correcting for this effect will result in a better match mY — K 
colour. 

It is also interesting to consider what these preliminary CB07 
models predict in terms of trends in the near-IR colours as a func- 
tion of age and metallicity. The Y — J colour exhibits similar be- 
haviour to optical colours, in that it increases with age. The J — H 
colour also increases slightly with age at mctalHcitics above OAZq, 
but it decreases with increasing age at lower metallicities. This ef- 
fect is even stronger for the H — K colour: a 0.2 Zq galaxy be- 
comes 0.2 magnitudes bluer as the galaxy ages from 5 Gyr to 15 
Gyr. The lower right panel of Fig. 5 shows the best combination 
of optical and near-IR colours for deriving age and metallicity de- 
duced from these preliminary CB07 models. As seen on this plot, 
the effect of increasing age is to move the galaxy along the g — r 
colour axis, without affecting Y — K colours. Conversely, increas- 
ing metallicity moves the galaxy along theY — K axis, with very 
little effect on the g — r colour. 



5.3.2 Non star-forming galaxies 

In Fig. 6 we compare the BC03 models with the colour distributions 
of galaxies in the non star-forming sample. Recall that all the mod- 
els on this plot have exponential declining star formation histories 
with r = 1 Gyr. The observed colours of non star-forming galax- 
ies have significantly less scatter than those of galaxies in the star- 
forming sample. We find that models with metallicities between 
OAZq andZQ reproduce the location of the observed galaxies in 
the the g — r versus u~ g colour-colour diagram, but the observed 
locus of the data in the J — H versus Y — J colour-colour plane 
cannot be understood with these models. The models also predict 
J — H colours that are too blue by about 0.1 magnitude. 

As seen from Fig. 7, the preliminary CB07 models predict that 
the near-IR colours of non-star-forming galaxies depend weakly 
on age at metallicities below 0.4 Zq . The age dependence is less 
pronounced than for the exponential declining models with r = 3 
Gyr that we used to model star-forming galaxies. This is due to the 
smaller fraction of luminous TP-AGB stars relative to RGB stars in 
the r = 1 Gyr models at the ages considered here. As was the case 
for the star-forming galaxies, we find that the preliminary version of 
the CB07 models do not match our data. These models are almost 
identical to the BC03 models at metallicities above 0.4Zq. The 
galaxies in our non star-forming sample are mainly massive, metal- 
rich systems, and the low-metallicity models are not applicable to 
them. 



5.4 Comparison with Maraston (2005) models 

5.4.1 Description of Maraston (2005) models 

One of the first papers to point out the influence of the short dura- 
tion thermally pulsating (TP-)AGB phase on the spectra of galaxies 
at A ~ was written by Maraston (1998). The most recent ver- 
sion of the Maraston models (Maraston 2005; hereafter M05) im- 
plements this phase by adopting a 'fuel consumption' approach. 
In Rettura et al. (2006) M05 model predictions were compared 
with those of BC03 and PEGASE.2 (Fioc & Rocca-Volmerange 
1997) for estimating the stellar mass content of early-t5?pe galaxies 
at « ~ 1. 

In the following sections we compare the predictions of CB07 
with those of M05. In Figs. 8 and 9, we plot the M05 models and 
the preliminary CB07 models on the top of our star-forming sam- 
ple and non star-forming sample in a series of optical and near- 
IR colour-colour diagrams. For the star-forming sample, the star 
formation history is an exponentially declining SFR with r = 3 
Gyr. The non star-forming sample is represented by an exponen- 
tially declining SFR with t = 1 Gyr. Results are shown at three 
different ages (5,10 and 15 Gyr), but only for half solar and so- 
lar metallicities (we do not have lower metalUcity predictions for 
M05). Figure 9 shows twice solar metallicity in addition. An ex- 
tinction of TV = 0.3 is already included in the models. The arrow 
on each panel indicates the amount of extinction corresponding to 
TV = 1 .4 as parametrised with a dust screen and the reddening curve 
of Cardelli et al. (1989), including the update for the near-UV given 
by O'Donnell (1994). 



5.4.2 Results from M05 models for the star-forming sample 

Figure 8 compares the predictions of CB07 with those of M05 for 
our star-forming sample, in the same way as described in the pre- 
vious section. As seen from this figure, both models agree with the 
data for the optical colours, M05 predicting slightly redder colours 
in g — r. In the J — H \s Y — J colour-colour plane, we notice 
that the preliminary CB07 models are considerably more sensitive 
to metalUcity. For both models, an extinction corresponding to tv 
= 1.4 is required to explain the reddest galaxies. This amount of 
extinction agrees with the measurement of the most dusty galax- 
ies in our star-forming sample. On the other hand, such amount is 
not necessary to explain the distribution of the optical colours. The 
predictions of M05 models match the data particularly well in the 
Y — K \s H — K colour-colour plane; in this plane, the prelim- 
inary CB07 models are offset m H — K. Both models exhibit a 
larger sensitivity to age at lower metallicity. The lower right panel 
of Fig. 8 shows that the M05 and CB07 models differ quite sig- 
nificantly in their predictions fotY — K colours. If one uses the 
preliminary CB07 model, one might infer that this colour is a good 
indicator of metallicity, but this is is not true for the M05 models. 
Based on this comparison, neither model is particularly favoured to 
explain the colour distribution of our star-forming sample. 



5.4.3 Results from M05 models for the non star-forming sample 

Figure 9 compares the predictions of CB07 with those of M05 for 
our non star-forming sample. We have included the 2Zq isometal- 
licity track as these galaxies have measured stellar metallities in 
the range 0.7^© to 1.5Zq. We note that M05 models require 
systematically higher metalUcities (between Zq and 2Zq) to fit 
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Figure 6. As in Fig. 4, except tliat the BC03 models are superimposed on our sample of non-star-forming galaxies. 



m 
< 



1.5 



1 .0 



0.5 



0.0 



: z 


= 0.02 Zs. 


.CB07 ; 




z 


= 0.2 Zs . . .. 






■ z 


= 0.4 Ze 




o 


z 


= ..^P 




cn 

<D 


: z 


= 2 Z, J* — ■ 




>^ 
X 






/ ■ 


1 



2 




u-g (AB) 



0.2 0.4 0.6 
Y-J (Vega) 




0.2 0.4 0.6 
H-K (Vega) 



< 




0.5 



0.0 



1.0 1.2 1.4 1.6 1.8 2.0 2.2 
Y-K (Vega) 



Figure 7. As in Fig. 4, except that the CB()7 models are superimposed on our sample of non-star-forming galaxies. 
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Figure 8. Colour-colour plots of M05 and CB07 composite stellar population models superimposed on our star-forming galaxies in black dots. The squares 
account for the M05 exponential star-forming with r = 3 Gyr models. The circles account for the CB07 exponential star-forming with t = 3 Gyr models. 
The grids are drawn for 3 metallicities and three ages (5, 10, 15 Gyrs): the colourful lines represent isometallicity colours of 0.5Zq, Zq. The amount of 
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Figure 9. As in Fig. 8, except that the M05 and preliminary CB07 models are superimposed on our sample of non-star-forming galaxies and that the metallicities 
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the observed colours than do the prehminary CB07 models (be- 
tween 0.5 Zq and Zq). The dust reddening should barely affect the 
colours as the non star-forming galaxies have a measured 2-band 
attenuation close to zero. Our conclusions regarding the fit of the 
optical and near-IR colours with both models are very similar to 
those reached for the star-forming sample. In optical colours, M05 
models predict too red g — r colours by approximately 0.2 magni- 
tudes. Both models predict too blue J — H colour by around 0. 1 
magnitude. The H ~ K colour is better reproduced by M05 mod- 
els. As pointed out in section 5.3.2, the age dependence is less pro- 
nounced than for the exponential declining models that are used to 
model the star-forming sample. The preliminary CB07 models pro- 
vide a better fit to the data mY — K colour; the M05 models requir- 
ing metallicities around 2Zq, which is larger than the metallicities 
estimated for the galaxies in our sample (Gallazzi et al. 2005). 
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6 IMPLICATIONS 

6.1 Can the age-metallicity degeneracy be broken? 

As discussed in Section 5.3, the g — r versus Y — K colour-colour 
plane may provide the ideal combination of colours to break the 
age-metallicity degeneracy. How well does this work in practice? 

In this section we analyze how galaxies with different mea- 
surements of stellar age and gas-phase metallicity populate the g— r 
versus Y — K colour-color plane. Fig. 10 shows the distribution of 
galaxies with different mean stellar ages. We have divided our sam- 
ple into four age classes, which contain equal numbers of galaxies. 
Blue points represent the youngest galaxies and red points the old- 
est. As seen from this Figure, the four age classes are distributed in 
horizontal slices, showing that as predicted by the models, the g — r 
colour is a good age indicator Figure 1 1 shows how galaxies with 
different gas-phase metallicities are distributed in the same colour- 
colour plane. Again, we have divided our sample in four classes: 
the blue crosses indicate the most metal-poor galaxies, while the 
red crosses are for the most metal rich ones. According to our pre- 
liminary version of the CB07 models, the Y — K colour should 
be very sensitive to the metallicity. The data reveals that galaxies 
in the lowest metallicity class do have predominantly blue Y — K 
colours, but galaxies in the medium to high metallicity classes all 
have similar Y — K colours. Dust is probably the main reason for 
this degeneracy. As seen from Fig. 5, reddening can influence the 
Y ~ K colours quite significantly, because the Y and K filters are 
distant in wavelength. 

In summary, our preliminary version of the CB07 models pre- 
dict that the location of galaxies in the g ~ r versus Y ~ K colour- 
colour plane allows one to estimate age and metallicity indepen- 
dently for a fixed amount of dust attenuation, but the data show that 
the method is not as clean as it might first appear. 



6.2 Estimating mass-to-light ratios 

When spectroscopic information is not available, the mass-to-light 
ratio of a galaxy can be estimated by comparing the predicted spec- 
tral energy distributions of model galaxies derived using stellar 
population synthesis codes to the observational data. The stellar 
mass of the galaxy is then obtained by multiplying its mass-to-light 
ratio by its observed luminosity. As discussed by De Jong (1996), 
it is important that these mass estimates be anchored by observa- 
tions in at least one photometric bandpass that is not too sensitive 



Figure 10. Distribution of galaxies with different mean stellar ages in the 
g — r versus Y — K colour-colour plane. The galaxies are divided in four 
equal classes according to their stellar age. Blue (triangles), green (filled 
circles), orange (stars) and red (crosses) indicate galaxies of increasing age. 
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Figure 11. As in Fig. 8, except galaxies are divided according to metallicity. 



to the presence of very young stars, or to metallicity and dust at- 
tenuation. The K-band is commonly regarded as the most reliable 
in this respect, because dust and young stars should have little ef- 
fect on the light emitted by a galaxy at wavelengths around 2 fim. 
Therefore, as pointed out by Maraston et al. (2006) and by Bruzual 
(2007), the influence of AGB stars on the model predictions for 
the AT-band luminosity have very important implications for the 
reliability of stellar mass estimates that make use of data in this 
wavelength range. In Fig. 12, we compare M,/Lk predicted by the 
BC03 and preliminary CB07 models as a function of the age of the 
galaxy. We show results for 5 different metallicities for a model 
with exponential declining SFR with r = 3 Gyr. As described in 
the previous section, this model provides a rather good fit to the 
optical colours of the star-forming galaxies in our sample. As can 
be seen from this Figure, at ages less than 1 Gyr, the mass-to-light 
ratio predicted by the CB07 models ranges from 30% (for 0.02^0 ) 
to 70% (for 2Zq) of the value predicted by the BC03 models. The 
difference between the two models decreases as a function of the 
age of the galaxy. At 10 Gyr, the typical age of the galaxies in our 
sample, the mass predicted by CB07 model ranges from 70% to 
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Figure 12. K-band stellar mass-to-light ratio of CB07 model galaxies di- 
vided by the K-band stellar mass-to-light ratio of BC03 model galaxies is 
plotted as a function of age for an exponential declining SFR and for 5 
metallicities from 0.02Zq to 2 Zq. 



100% of the value predicted by the BC03 model. We conclude that 
even for present-day star-forming galaxies, there are significant un- 
certainties on stellar masses derived from A'-band luminosities that 
depend on how AGB stars are treated in the population synthesis 
models. The stellar masses predicted by the CB07 models are con- 
siderably smaller that those predicted by the BC03 models, the ef- 
fect being more important at high redshifts and for low metallicity 
galaxies. 



7 SUMMARY 

We have used a combination of empirical and theoretical tech- 
niques to interpret the near-IR colours of 5800 galaxies drawn from 
the SDSS main spectroscopic sample. This study focuses on the in- 
ner regions of galaxies that are sampled by the 3 arcsecond SDSS 
fibre spectra. In the first part of our analysis, we study correlations 
between near-IR colours and physical parameters derived from the 
spectra, which include the specific SFR, the stellar age, the metal- 
licity, the dust attenuation and the axis ratio of the galaxy. All cor- 
relations are analyzed using galaxy samples that are matched in 
redshift, stellar mass and concentration parameter. We analyze star- 
forming and non-star-forming galaxies separately. In the second 
part of our analysis, we compare the near-IR colours of the galaxies 
in our sample to the predictions of stellar population models. Our 
main conclusions are as follows: 

• Whereas more strongly star-forming galaxies have bluer op- 
tical colours, the opposite is true at near-IR wavelengths: galaxies 
with higher values of SFR/ A/, have redder near-IR colours. This 
result agrees well with the predictions of models in which TP- AGB 
stars dominate the H and i^-band light of a galaxy following a 
burst of star formation. 

• We find a surprisingly strong correlation between the near-IR 
colours and dust attenuation measured from the B aimer decrement. 
However, the near-IR colours do not correlate with the axis ratios of 
the galaxies. This suggests that the correlation with dust attenuation 
arises because TP-AGB stars are the main source of dust in the 
galaxy. 

• We compare the near-IR colours of the galaxies in our sample 



with the colours predicted by the BC03 code, a preliminary version 
of the CB07 code and the M05 code. The preliminary CB07 code 
provides better qualitative agreement with the data than BC03, in 
that it predicts that more strongly star-forming galaxies have redder 
near-IR colours. However, this effect occurs only at low metallic- 
ities. The observed spread in the the near-IR colours of the star- 
forming galaxies is still unexplained by the current models. The 
M05 models agree well with the data in the Y — K \s H ~ K 
colour-colour plane, but do not provide as a good a match for other 
colour combination. 

• The preliminary version of the CB07 model suggests that both 
the age and the metallicity of a galaxy can be estimated from its 
location in the g ~ r versus Y — K colour-colour plane. We have 
tested this using our data and we find that although age can be quite 
accurately estimated from the g — r colour, the Y — K colour can 
only be used to distinguish the lowest metallicity galaxies. 

• Even for present-day star-forming galaxies, the if-band mass- 
to-light ratios predicted by the CB07 models can be significantly 
(up to 30%) smaller than those predicted by the BC03 models. 
This systematic uncertainty in the mass estimates is comparable 
to the differences obtained when stellar masses are estimated using 
a variety of different methods, for example fitting stellar population 
synthesis models to broadband colors or analysis of spectroscopic 
features (Drory et al. 2004). 

Throughout this paper, we have been careful to stress that the 
model comparisons have been made using a preliminary version 
of the CB07 model. The published version of the model should 
include an improved library of spectra of TP-AGB stars, as well 
as new stellar evolutionary tracks for the pre-AGB phase. This is 
likely to improve the agreement with the data. The hope is that 
the sample of nearby galaxies with high quality photometric and 
spectroscopic data from UKIDSS and SDSS presented in this paper 
will continue to serve as a means of testing and calibrating future 
stellar population synthesis models. 
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